Analytical description of a trapped semi-ideal Bose gas at finite temperature 
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Present experiments with Bose condensed gases can be largely described by a semi-ideal two-gas 
model. In this model, the condensate is influenced only by the mean-field repulsion among condensed 
atoms, while the thermal cloud is considered an ideal gas confined by an effective potential that 
consists of the external trap and the mean-field repulsion by the condensate. This simple, intuitive 
model provides explicit analytical expressions for the density distributions of the condensate and 
the thermal component. It describes the reduction of the condensate fraction relative to that of an 
ideal gas as a consequence of the positive chemical potential due to interactions in the condensate. 
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Recent experiments with Bose-Einstein condensed al- 
kali gases (cf. e.g. have been successfully described 
by different approximation schemes. The diluteness of 
these gases, expressed as na 3 <C 1 where a is the two- 
body scattering length, allows a description in terms of an 
effective mean-field theory. At zero temperature the con- 
densate can be described by the Gross-Pitaevskii equa- 
tion a nonlinear Schrodinger equation. In present 
experiments interactions are strong enough so that the 
chemical potential /i is much larger than the level spac- 
ings Tiu) x , hujy, and tuv z of the external harmonic poten- 
tial. Therefore the use of the Thomas-Fermi approxima- 
tion which neglects the kinetic energy of the con- 
densate, has lead to a tremendous simplification in the 
understanding of these experiments. 

Recently, much theoretical and experimental atten- 
tion has been focused on finite temperature properties 
of Bose-condensed gases. A semiclassical Hartree-Fock- 
Popov description has produced excellent agreement with 
experimentally measured equilibrium quantities [Q . Fur- 
ther, unlike in a spatially homogeneous gas, it was found 
that collective elementary excitations contribute very lit- 
tle to thermodynamic properties due to the relatively 
small spatial overlap between condensate and thermal 
component in a harmonic trap. Thus, even a semiclas- 
sical Hartree-Fock approximation is a fairly accurate de- 
scription of these gases, except for a small temperature 
range close to the critical temperature. For example, ex- 
cellent agreement has been achieved by a comparison of 
Hartree-Fock density profiles with the results of a com- 
prehensive Path Integral Monte Carlo calculation || . 

In this paper we introduce a further simplification 
of the self-consistent Hartree-Fock model by neglecting 
atomic interactions in the thermal component of the gas. 
Like the introduction of the Thomas-Fermi approxima- 
tion Q the reduction to Eqs. (^|)-@ and ([lj) has the 
important advantage of providing an explicit analytical 
description of a trapped Bose gas at finite temperature. 



We consider atoms that are trapped in a generally 
anisotropic harmonic potential. Under the appropriate 
rescaling of energy and lengths, the harmonic potential 
has the form V(r) = r 2 /2. Here, the dimensionless carte- 



sian coordinates of r are given in units of ^Jhuj/mu}^, y z , 

with u; = {uJxUjyUJz) 1 / 3 . Other lengths, like the thermal 
wavelength At and the scattering length a, are scaled in 
units of the natural length yjfi/muj of a harmonic oscil- 
lator with angular frequency lu. The energy V(r) as well 
as other energies are implicitly given in units of Tilo. 

Then, in the self-consistent Hartree-Fock model, the 
densities of the condensate and the thermal component 
are given as solutions of the two coupled equations PJ7l 



no(r) = " ~ " 2UnT{r) 0(» - r 2 /2 - 2Un T (r)) (1) 

nT ( r ) = _L 53/2 ( e -(r 2 /2+2C/ [„o(r)+n T (r)]- M )/fc B T^ _ (2) 

Here, the Thomas-Fermi approximation has been applied 
which neglects the kinetic energy of the condensate and 
implicitly finite size effects. The strength of the atomic 
interactions is given by U = Ana. We use the thermal 
wavelength A^ = -\/27r/fcsT and the Bose function g 3 / 2 is 
defined by g a (z) — JZj—i z3 1 J a ■ T ne chemical potential 
fj, is determined by the constraint of a given total atom 
number 



N= dr [no(r)+nr(r)]. 



(3) 



For further reference, Eqs. (|l|) - (||) will be denoted as 
the interacting model. Even though it represents already 
a relatively concise description, these equations still have 
to be solved self-consistently. This is done numerically 
by an iterative procedure Q . 

If one further neglects the mean-field repulsion from 
non-condensed atoms Eqs. ([j]) and (|^) are solved by the 
explicit relations 
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with the chemical potential fi given by 
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These relations summarize a simple model of the partly 
condensed gas which we denote as the semi-ideal model. 
The condensate density uq is that of a zero temperature 
condensate with Nq atoms, while the density n-r of the 
non-condensed cloud is that of an ideal gas of bosons 
confined in the combination of the external potential and 
the repulsive mean-field potential due to the condensate 
atoms. Its wide range of validity has been confirmed by 
a numerical analysis in Ref. ||. However, in order to 
make it an intuitively appealing, usable description of 
a trapped Bose gas, the condensate fraction has to be 
determined analytically as a function of temperature. 

For later convenience let us introduce the reduced 
chemical potential /2 
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Here, the critical temperature T c is that of a trapped 
ideal gas @,0 k B T c = [iV/C(3)] 1/3 where g a (l) = ((a) 
was expressed in terms of the Riemann Zeta function. 
The reduced chemical potential depends on the scaling 
parameter 77 H 
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which describes the strength of the atomic interactions 
within the condensate. The scaling parameter is indepen- 
dent of the system size when the thermodynamic limit is 
taken in the usual way (N — > 00, ui — > 0, Nlu 3 — const.). 
Due to its relatively weak TV 1 / 6 dependence, rj assumes a 
value close to 0.3 in most recent experiments 0. 

The condensate fraction is determined by integration 
over the thermal density distribution. For T < T c , the re- 
sult can be written in terms of incomplete Gamma func- 
tions @ 
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Eq. (|^) is still an implicit expression since it depends 
again on the reduced chemical potential. Therefore, R 
and Nq/N are given as simultaneous solutions of Eqs. (0) 
and (||) . For temperatures T > T c a condensate does not 
exist and the reduced chemical potential fl < is trivially 
defined by the single equation g 3 (e») = ((3)(T/T C )- 3 . 

Let us consider the results of the semi-ideal model in 
the limits of small and large /2, corresponding to T « T c 
and T « 0. Eq. (ft) is therefore rewritten 
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in terms of the function 
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(11) 



where the reduced energy e = e/ksT has been intro- 
duced. The function f(jp) can be expressed as a power 
series in y/Jp,, where ft runs from (T = T c ) to 00 
(T = 0). Truncating the series after the lowest nontrivial 
order yields 
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The first two terms of Eq. ( |l2|) correspond to the conden- 
sate fraction of an ideal gas, whereas the third describes 
the influence of the condensate repulsion to lowest order 
in /2. This expansion adequately describes the solution of 
Eq. ([)]) over the entire range of temperatures, provided 
that 77 <C 1. A further simplification is derived by solving 
Eq. (O) to lowest order in 77, arriving at 
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Within Eqs. (|12|) and (|13|), the only effect of the con- 
densate interactions is to shift the chemical potential to 
positive values, raising the bottom of the effective poten- 
tial which the thermal atoms experience. The effect of 
changing the shape of the potential enters only at higher 
orders of the expansion. Eq. ( |l2| ) is thus a good ap- 
proximation of Eq. (^J) if most of the thermal atoms are 
spatially separated from the condensate, a condition sat- 
isfied for temperatures close to T c and for small 77. 

Eq. (12) is not expected to be a good description for 
extremely low temperatures. The understanding of this 
limit is facilitated by f(x) obeying the differential equa- 
tion 



d^ f - f - 2 ^- 



(14) 
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The function f(x) has to grow slower than x 3 / 2 for large 
values of x. Otherwise, the power series in Eq. (|io| ) would 
not converge. As a consequence, the second derivative in 
Eq. (^) vanishes for large x. Thus, the asymptotic limit 
of the condensate fraction for low temperatures is given 
by f(x) = 2y/x or 
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However, it has been shown [Q] that the leading term of 
the uncondensed fraction around zero temperature scales 
like 
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due to quasi-particle contributions which are included in 
neither the semi-ideal nor the interacting Hartree-Fock 
treatment. 

According to Fig. |l|, where ?y = 0.31 has been assumed, 
the high temperature expansions Eqs. ( |l2[ ) and to some 
extent Eq. (Jl^) are very good approximations of the exact 
condensate fraction Eq. (0) of the semi-ideal model over 
the whole range of temperatures T < T c . In contrast, 
we have found that the validity of the low temperature 
limits of Eqs. ( |l5| ) and ( |l6| ) are restricted to a very small 
range around zero temperature which makes them nearly 
irrelevant for practical purposes. 

The condensate fraction given by the approximate so- 
lution of the semi-ideal model in Eq. ( |l2|) is compared 
in Fig. p| with the result of an iterative solution of the 
interacting model (Eqs. (|l|)-(§)). The difference between 
the curves is small except near the critical temperature, 
where the interactions among thermal atoms shift the 
critical temperature jl3|]. This shift is not exhibited by 
the semi-ideal model. The wide range of validity of the 
semi-ideal model has also been confirmed by a numer- 
ical comparison with a Hartree-Fock-Popov calculation 
The latter includes collective excitations, which are 
ignored in the models shown in Fig. |^. These lead to 
an additional but negligible decrease of the condensate 
fraction [Q. 

A direct comparison between the density distributions 
of the different models is given in Fig. || for T/T c = 0.5. 
Here, the Thomas-Fermi solution of the condensate den- 
sity of Eq. (Q) has been replaced by a solution of the full 
Gross-Pitaevskii equation. The main implication of the 
included kinetic energy of the condensate is a smoothed 
out condensate surface. Within the Thomas-Fermi ap- 
proximation almost no difference is seen between the two 
curves. Even though the unaccounted existence of collec- 
tive excitations changes the densities of condensate and 
thermal component close to the center of the trap by a 
significant amount, their influence on the total density is 
negligibly small for most temperatures Q . 

In recent experiments, trapped Bose gases have been 
probed in-situ using either non-destructive dispersive 



ima ging techniques |T^] , or off-resonant absorption imag- 
ing |p.5| . These methods are important new tools in the 
understanding of trapped Bose gases, with distinct ad- 
vantages over the previously used time of flight measure- 
ments of expanding atom clouds [jiJU . We therefore ex- 
pect that the results of this paper will become useful 
for future in-situ experiments with trapped Bose gases. 
These new optical detection techniques measure column 
or line density profiles; we therefore determine these pro- 
files from Eqs. (|J) and (|^) by integration over one or two 
dimensions. 

In cylindrical coordinates, where the probe beam is 
directed along the z-axis, the according column densities 
are given by 
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where /i(p) = (fi — /9 2 /2)/fceT. Equivalently, we obtain 
the line densities 

no(z) = 7r ^-^ /2)2 6(»-z 2 /2) (19) 



X 5 T l2C(5/2)- 55/2 (e-PW) fi(z)>0 

with p,(z) = (fi- z 2 /2)/k B T. 

We therefore conclude that the proposed semi-ideal 
model of Eqs. @-@ and ( |l2| ) represents a rather good 
description of a trapped Bose gas at finite temperature 
provided that r] is not larger than the values r\ ~ 0.3 
to 0.4 achieved in recent experiments. Its mathematical 
simplicity gives a much clearer conceptual picture than 
the self-consistent approaches that have been used be- 
fore. The observed agreement of the semi-ideal model 
with a self-consistent Hartree-Fock calculation has con- 
firmed that the density of the thermal component is too 
low to necessitate an interacting gas description. Instead, 
the observed depletion of the condensate fraction arises 
from condensate interactions, leading to a positive value 
of the chemical potential or equivalently to a raised bot- 
tom of the effective potential in which the thermal com- 
ponent resides. It is well described by Eq. (|l2|), which is a 
low order expansion in p,. In contrast, the changed shape 
of the effective potential does not influence the thermal 
population in a noticeable way. 
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FIG. 1. Condensate fraction as a function of temperature, 
given by the semi-ideal model. The dashed line shows the 
exact solution Eq. (^|) of the semi-ideal model. The solid line 
represents its high-temperature limit Eq. ([l2]). The explicit 
expression of Eq. ( |l3[ ) is given by the dotted line. This plot 
demonstrates the validity of Eq. (jl^) over the whole range of 
temperatures T <T C . The explicit approximation of Eq. ( p^ ) 
fails only in the vicinity of the critical temperature (shown in 
the inset). 
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FIG. 2. Condensate fraction plotted as a function of tem- 
perature, comparing the discussed models. The ideal gas 
condensate fraction is given by the dashed-dotted line. The 
dotted line is a numerical solution of the Hartree-Fock Eqs. 
(S) - ©- The solid line is obtained from the analytical ap- 
proximation Eq. (^), corresponding to the solid line of Fig. 
hi Except near the critical temperature (shown in the inset), 
the semi-ideal model well approximates the solution of the 
interacting model. The critical temperature T c refers to the 
case of an ideal gas. 




FIG. 3. A comparison of the density distributions at 
T/T c — 0.5 predicted by the two discussed models. The up- 
per solid line shows the density of the thermal component 
as given by Eq. (^). In comparison, the lower solid line has 
been calculated numerically using Eq. (^). However, a nu- 
merical solution of the full Gross-Pitaevskii equation has been 
used instead of the Thomas- Fermi approximation of Eq. (|l|) . 
A small difference exists at the edge of the condensate due 
to the neglect of the kinetic energy in the semi-ideal model. 
However, the total density is only weakly affected by this ef- 
fect, as is seen by the two dashed-dotted lines. Again, the 
upper dashed-dotted line refers to the total density given by 
the semi-ideal model while the lower curve corresponds to 
the self-consistent Hartree-Fock approximation. Here, a total 
atom number of N = 5 x 10 6 was assumed. 
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